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Noise and Delay Induced on Dynamics of the 
Aeroelasticity Model 

Zaitang Huang 


Abstract. The purpose of this paper is to investigate the 
stochastic bifurcation and stability problem of the Aeroelasticit 
-y of two-dimensional supersonic lifting surfaces with delay 
term. Applying Hopf bifurcation theory, Lyapunov exponent 
and invariant measure theory, we analyze the D- and 
P-bifurcation of the stochastic system. The analysis is based on 
the reduction of the infinite-dimensional problem to one 
described on a two-dimensional stochastic center manifold. 
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dynamical system provides a very powerful mathematical 
tool for understanding the limiting behavior of stochastic 
system. It has been applied to engineering, respiratory 
physiology, chemical plants, mechanical systems, fluid 
dynamics, secure communications, economics and biological 
systems[16-21]. Our purpose in this paper is to investigate 
the stochastic dynamical behavior for the system (2.1) by 
applying the singular boundary value theory, Lyapunov 
exponent and the invariant measure theory, the direction of 
the Hopf bifurcation and the stability of bifurcating periodic 
solutions are also determined. 


I. Introduction 

Because of its evident practical importance, the study of 
the flutter instability of flight vehicle constitutes an essential 
prerequisite in their design process. The flutter instability can 
jeopardize aircraft performance and dramatically affect its 
survivability. Moreover, the tendency of increasing structural 
flexibility and maximum operating speed increases the 
likelihood of the flutter occurrence within the aircraft 
operational envelope. As a result of the considerable 
importance of this problem, a great deal of research activity 
devoted to the aeroelastic active control and flutter 
suppression of flight vehicles was carried out. In this sense, 
the reader is referred to a sequence of issues in Refs[l-5], 
where valuable contributions to this topic have been supplied. 

As it clearly appears, within this problem, two principal 
issues deserve special attention: 1) increase, without weight 
penalties, of the flutter speed, and 2) possibilities to convert 
unstable limit cycles into stable ones. While the achievement 
of 1) can result in the expansion of the flight envelope, the 
conversions mentioned in 2) would make it possible to 
operate in close proximity of the flutter boundary without the 
danger of encountering the catastrophic flutter instability, but 
in the worst possible scenario, crossing the flutter boundary 
that features a benign character. In contrast to the 
catastrophic flutter boundary in which case the amplitude of 
oscillations increases exponentially, in the case of benign 
flutter boundary, monotonic increase of the oscillation 
amplitude occurs in cases 1) and 2) respectively. And, as a 
result, the failure can occur only by fatigue. It clearly appears 
that both issues 1) and 2) are related to controlling Hopf 
bifurcations. In particular, issue 1) implies increase of the 
stability of an equilibrium and noise of the occurrence of 
Hopf bifurcations[6-9]whereas issue 2) is related to 
controlling Hopf bifurcations once a periodic vibration has 
been initiated [10- 15]. Recently, the theory of random 
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II. Stochastic aeroelastic model 


This investigation is based on a geometrical and 
aerodynamic nonlinear model of a wing section of the 
high-speed aircraft incorporating an active control capability. 
As concerns the nonlinear unsteady aerodynamic lift and 
moment, these are obtained through the integration of the 
pressure difference and of its moment with respect to the 
pitching axis, respectively, on the upper and lower surfaces of 
the airfoil. To this end, the third-order approximation of the 
piston theory aerodynamics[10-15](PTA), as given by 

p{x, t) = p yj {l + k(v 7 I a Jr + [k{k + 1) / 4] [(y . / a x )/] 2 

Hi 


+[K(K+l)/l2][(vJaJrT 

is considered, where K is polytropic gas coefficient. Here in 


dw 


dw 


V z=- 


— + U„ — 


dt 


dx 


5gn(z) 


denotes the downwash velocity normal to the lifting 
surface a 2 = KJ) v / p x where sgn(z) assumes the value 1 or 


— 1 for z > 0 and z < 0 , respectively. In addition, 
w(t) = h(t) + a(t)(x-bx 0 ) 

denotes the transversal displacement of the elastic 
surface; X () (= V ea ) is the dimensionless streamwise position 


of the pitch axis measured from the leading edge; [). r , p x , 

U x and Cl y are the pressure, the air density, the airflow speed, 
and the speed of sound of the undisturbed flow, respectively; 
and y — M y / M y _ — 1 is an aerodynamic correction factor 

that enables one to extend the validity of the PTA to the entire 
low-supersonic/hypersonic-speed range. 

As there also exist many stochastic factors affecting and 
disturbing the realistic environment considering the change 
of the twist angle about the pitch axis. We think it is 
reasonable and necessary to add random terms in the 
aeroelastic model. In the context of the inclusion of the 
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structural and aerodynamic nonlinearities, of the linear and 
nonlinear controls and of the associated noise and time delay, 
in conjunction with the typical cross section with pitch-and 
plunge degrees of freedom, the dimensionless stochastic 
aeroelastic equations representing an extension of those in 
Refs [18,26, 28] are written as 


x + % a a + 2£ h 


yV J 


V V J 


x = L{t), 


^-x(t) + d + — T a + — -Ba 3 =M(t)-^-La(t-r) 








VJ/ 

- — f a 3 (t - t) + oj (i + d)%{t) + (7 2 Tj{t ), 


where 


( 2 . 1 ) 


Lit) = ] - — {12a + 5 A Mli\ + K)y 2 a 3 

12 juM x 1 A 

+ 12 [x + aib-v ea )/b]}. 


M ^— T r Z r r {l^-yJa + S,M 2 Jb-vJ 
12 MM X 

(1 + K)y 3 a 2 + 4[3 ib - v ea )x + d(4Zr - 6b 

K.+y. )/6]}. 

and xit) = hit) / bill), CC it) is the twist angle about the 
pitch axis, ^it) is the multiplicative random excitation and 
r/it) is the external random excitation directly (namely 
additive random). d,if) and r/it) are independent, in 
possession of zero mean value and standard variance Gauss 
white noises, i.e. E[%{t)\ = E\r/it)\ = 0, £’[^(f)^(f + r)] 

= Sir), E[r/it)r/it + r)] = Sir), E[dit)r/it + r)] = 0. An- 
d (oj , a 2 ) is the intensities of the white noise, and Lit) and 
M it) denote the dimensionless aerodynamic lift and 
moment, respectively. The meaning of the remaining 
parameters can be found in the nomenclature (see also 
Refs[l-5, 10-15]). In Eq. (2.1), the parameter B identifies the 
nature of the structural nonlinearity of the system in the sense 
that, corresponding to B <0 or B > 0 , the structural 
nonlinearities are soft or hard, respectively, whereas for 
6 = 0 the system is structurally linear. The linear and nonli- 
near active controls are given in terms of two normalized 
control gain parameters and V F, , respectively. 

A mathematical model is generally the first approximation 
of the considered real system. More realistic models should 
include some of the stochastic factors affecting and past 
states of the system, that is, the model should include noise 
and time delay. The noise and time delay in control can occur 
either beyond our will or it can be designed as to modify the 
performance of the system. For this reason, as a necessary 
prerequisite, a good understanding of its effects on the flutter 
instability boundary and its character (benign or catastrophic) 
is required. 

To capture the effect of noise %(t) , r/it) and time delay 
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T , introduced in the related terms and ® 2 , letx = x x , 

a = X 2 ,X = X 3 ,a = X 4 ,X 2t = X 2 it — r). Then, one can 

rewrite Eqs.(2.1) as a set four first-order differential equatio- 
ns: 

i, = x 3 , 

x 2 (t) = x 4 , 

x 3 it) = a x x j + a 2 x 2 + a 3 x 3 + a 4 x 4 + a 3 x 3 2 + a 6 x 2t 

< 

+ a 7 x 3 , + ie x x 3 + e 2 x 4 )<^it) + e 3 r/it), 
x 4 it) =b l x l +b 2 x 2 +b 3 x 3 +b 4 x 4 +b 3 x 3 2 +b 6 x 2t +b 1 x i 2t 
+ id x x 3 + d 2 x 4 )^it) + d 3 r/it), 

(2.2) 

where a 6 = a 6c + a 6 is the bifurcation parameter, all of the 

coefficients that are provided in[12]. It is obvious that there 
exist a unique equilibrium point Q(0, 0, 0, 0). 

For convenience in the following analysis, rewrite Eqs.(2.1) 
in the vector form: 

xit) = Axit) + Bxit — t) + Fixit), xit — t), d, it ), r/it)), 
where X, F e jR 4 , A and 6 are 4 x 4 matrices. A, B, and 
F are given by 



"0 

0 

1 

o' 


"0 

0 

0 

0" 


0 

0 
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0 

0 

0 

0 

A = 

a { 


a 3 

a 4 

B = 

0 

0 

a 6 

0 


A 

b i 

h 

b 4 _ 


0 

K 

0 

0 


0 

0 

F = 

a 5 x 2 +a 7 x 2 it-r) + ie x x 3 + e 2 x 4 )^(t) + e 3 r/it) 
b 3 x\ +b 7 x\it -r) + id x x 3 + d 2 x 4 )J;it) + d 3 r/it) 
respectively. 

Hopf bifurcation has been extensively studied using many 
different methods[23-28] for example, Lyapunov quantity 
used in the context of the supersonic panel flutter[13-15] 
where the effects of structural, aerodynamical, and physical 
nonlinearities have been incorporated. In Refs[ll,15] the 
dynamic behavior of the system without noise and time delay 
in the control was studied in the vicinity of a Hopf bifurcation 
critical point. In particular, the effect of the active control on 
the character of the flutter boundary (where the Jacobian has 
a purely imaginary pair) is investigated. It is shown that for 
different flight speeds, stable (unstable) equilibrium and 
stable (unstable) limit cycles exist. 

The effect of the noise and time delay involved in the 
feedback control will be considered in this paper. Nonlinear 
systems involving time delay have been studied by many 
authors [5, 12, 13]. In the past two decades, there has been 
rapidly growing interest in bifurcation control[l-5, 10-15] 
There are a wide variety of promising potential applications 
of bifurcation and chaos control. In general, the aim of 
bifurcation control is to design a controller such that the 
bifurcation characteristics of a nonlinear system undergoing 
bifurcations can be modified to achieve some desirable 
dynamical behaviors, such as changing a subcritical Hopf 
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bifurcation to supercritical, eliminating chaotic motions, etc. 
In this context, many applications have been found, for 
example, in the areas of mechanical systems, fluid dynamics, 
biological systems, and secure communications. Although 
effects of delay time on the Aeroelastic model have been 
extensively investigated [1-5,10-15], there have been no such 
studies on the effect of multiplicative noise, to the authors’ 
knowledge’ . We are interested in the stochastic bifurcation, 
which is one of the interesting phenomena induced by noise 
(Refs. [16-21,25-29], related references therein). The main 
attention is focused on Hopf bifurcation. 

As the first step, we analyze the stability of the trivial 
solution of the linearized system of Eq.(2.3), which is given 
by 

x(t) = Ax(t) + Bx(t -r),re IK. ( 2 . 4 ) 

The characteristic function can be obtained by substituting 
the trial solution x(t) = ce/' , where cis a constant vector, 
into the linear part to find 

A(A,a 6 ) = det(AI-A-Be Ar ) = 0, (2.5) 

where I denotes the identity matrix. It can be shown[5] that 
the number of the eigenvalues of the characteristic equation 
(2.5) with negative real parts, counting multiplicities, can ch- 
ange only when the eigenvalues become pure imaginary pairs 
as the time delay T and the components of A and B are 
varied. 

It is seen from Eq.(2.5) that when 5a, (b 2 + b 6 )^ b x (a~, + 


a 6 ) none of the roots of A(A,Cl () ) is zero. Thus, the trivial 


equilibrium X = 0 becomes unstable only when Eq.(2.5) has 
at least a pair of purely imaginary roots A — +0)i ( i is the 
imaginary unit), at which a Hopf bifurcation occurs. 

To obtain the explicit analytical expressions for the 
stability condition of Hopf bifurcation solutions, system (2.1) 
should be reduced to its center manifold[12, 30-32]. While 
studying 

the critical infinite dimensional problem on a two-dimension- 
al stochastic center manifold, we express the delay stochastic 
equation as an abstract stochastic evolution equation on 
complete probability space. By the centre manifold theorem 
and Hopf bifurcations [12,30-32], we obtain the equation (2.1) 
of the stochastic center manifold: 


y(t) 


0 

-CO 


y + H(y), x(t - t) = O(-r)y(0, ^ 

(2.6) 


x(t) = G>(0)y(t), 

where H(y) represents the nonlinear terms contributed from 
the original system to the stochastic center manifold. 

The lowest-order nonlinear terms of the stochastic center 
manifold, needed to determine the solutions, are 

H 3 (y) = i//(0)F(Oy) = ® T (0) 


0 

0 

fl 5 (0(0)y) 3 2 +a 7 (0(-r)y) 3 2 +( ei 0(0)y + e2 0(0) 
y)%(t) + e 3 r/(t) 

b 5 m0)y)l+b i mT)y) 3 2 + (d l O(0)y + d 2 O(0) 
_y)4(t)+d 3 Tj(t) 

c nyi + c i 2 yfy 2 + c !3>'i 3 ; 2 2 +Ci 4 y 2 3 +(eiyi+e 2 y 2 ) 

£(t) + e 3 77(f) 

^ 21^1 ~^ C 22^1 3^2 ^^ 23 ^ 1^2 ^^ 24^2 (^ 1 3 ^ 1 ”^ ^ 2 ^ 2 ^ 

_£(t) +d 3 T](t) 
where 
O>(0) = 


COS 0)0 

sin 0)0 


L, cos 0)0 + coL 2 sin co0 

L, sin o)0 - o)L 7 cos o)0 

h 

A, 


-a) sin0 

0) COS 0)0 


0 )(a)L 2 cos o)0 - Lj sin o>9) o)( o)L 2 cos co0 + L x sin co0) 

L A, 

L 0 





L, cos o)0 + L 4 sin o)0 

sin 0)0 - L 4 cos 0)0 


M 

M 


L 5 cos o)0 + L 6 sin o)0 

L 5 sin o)0 - L 6 cos co0 


M 

M 

5 

L 7 cos o)0 + l ^ sin o)0 

Lj sin 0)0 - 1^ cos 0)0 


M 

M 


N l cos 0)0 - N., sin o)0 

A, sin o)0 + N , cos o0 


where the explicit expressions of L t (i = 1,2,.., 8) and N 
are also provided in[12], and N } and N 7 can be obtained 

from the relation ( l l i , A> ) 

= I , expressed in terms of 0), T 

and the coefficients a n 

bj,d n and e, in Eqs.(2.1). The 


lengthy expressions of C jf , e t and d t are omitted here. 

Therefore we obtain the equation (2.1) of the stochastic 
center manifold: 

A (0 = +c io> ; i + c ioo > ; 2 + c n>'i 3 +c i 2 y 2 y 2 + 

c^yjl +c u yl + (e 1 y 1 +e 2 y 2 )Z(t) + e 3 Ji(t), 

< 

y 2 ( 0 = -ay, + c 20 yi + c 200 y 2 + c 2 i y i 3 + c 22 y\ y 2 + 

c 22 y{y\ +c 2A y\+(d x y x +d 2 y 2 )^(t) + d 3 ri(t). 
We set the coordinate transformation Vj = r COS 0, y 4 = 
r sin 9, and by substituting the variable in (2.7), we obtain 
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r{t) = rc l0 cos 2 9 + r(c m + c 20 ) sin 9 cos 9 + rc 200 
sin 2 0r 3 [c n cos 4 9 + {c n +c 21 )cos 3 9sin9 
+ (c 13 + c 22 )cos 2 (9sin 2 9 + {c u + c 23 ) cost? 
sin 3 9 + c 24 sin 4 9 ] + r[e x cos 2 9 + (e 2 + d x ) 

cos 9 sin 9 + d 2 sin" 9]g(t) + [e 3 cos 9 + d 3 
sin 9]ij(t), 

i 9{t) = -co + c 20 cos 2 9 + (c 200 -c 10 ) sin 9 cos 9 
— c 100 sin 2 9 + r 2 [c 2l cos 4 9 + (c 21 -c 12 ) 
cos 3 9 sin 9 + (c 23 + c 12 ) cos 2 9 sin 2 9 + 

(c 24 +c 213 )cos(9sin 3 9 + c 14 sin 4 0] + [j, 
cos" 9 + (d 2 -e x )cos9sin9-e 2 sin" 9 ] 

%{t) + -[J 3 cos 9 + <? 3 sin 9]rj(t). 


(2.8) 

It is difficult to calculate the exact solution for system (2.8) 
today. According to the Khasminskii limit theorem, when the 
intensities of the white noises (e j , d i )(i = 1,2,3) is small en- 


enough, the response process {r(t), 9(t)} weakly converged 


to the two-dimensional Markov diffusion process [26-29]. 
Through the stochastic averaging method, stochastic 
differential equations (2.9) are obtained 


\dr = m r dt + cr, ] dW r + a 12 dW g , 
\cl9 = m g dt + a 21 dW r + a 22 dW g , 


(2.9) 


where W r (t ) and W g are the independent and standard 

Wiener processes. As for the twodimensional diffusion 
process, it is necessary to calculate its two-dimensional 
transition probability density. There is no general and right 
method for the calculation. As for the concrete system, we 
could finish the calculation with some special techniques. 

Set the parameters as follows: 

M ~ “ (^to "*"^200)’ Mo — ~tt> + — (^20 +^ioo)> 

jU 0 = 5s x + 5d 2 + 3d[ + 3c 2 t foc 2 d x — 2 e x d 0 , 

ju 3 = — (c 3 + d 3 ), jU 4 = 3e~ + 3d-, + + d x 

+ 2 e 2 d x + 2e x d 2 , 

jU 5 = — + d 0 )(d[ — c 2 ), /fg = + d + 3e 2 

+ 3d 2 - 2e 2 d x - 2 e x d 2 , 

M — 3c u +3 c 21 +c 13 +c 22 ,/f 8 — 3 c 21 —3c h + c 23 —c 12 . 


Under the condition cr 12 = <j\ A 0, we rewrote system (2. 
9) as follows 


dr = 


M + 


M 


r + — + — r 3 


dt + 


ju 3 + — r 2 

8 


dW r + (rju 5 ) 2 dW g , 

\ 


d9 


Mr 2 

Mo + yr" 


v 

dW a . 


dt + (r /u 5 )~ dW r + 


( hL 

2 + 


V ' 


8 


( 2 . 10 ) 

From the diffusion matrix, we can find that the averages 
amplified r(t) is a one-dimensional Markov diffusion 

process when cr 12 = <r 21 = 0, i.e. e x +d 2 = 0, or d x ~e 2 = 0 
.Thus we have the equation as following 

7 ■■ N \ 


dr = 


M + 


M 

8 


r + — + — r 3 


8 


dt 


( 2 . 11 ) 


+ 


ju 3 + — r 2 

3 8 


dW,. 


This is an efficient method to obtain the critical point of 
stochastic bifurcation through analyzing the change of 
stability of the averaging amplitude r(t ) in the meaning of 
probability. 


III. Stochastic D-bifurcation 

In the section. We will see how the introduction of 
randomness change the stochastic behavior significantly 
from both the dynamical and phenomenological points of 
view[26-29]. 

Theorem 3.1 (D-Bifurcation) When //. = 0, jU 7 = 0. 
Then the delayed stochastic system (2.2) undergoes a 
D-bifurcation, at the parameter value 16//, +2 fl, = /J 4 . 
But the stochastic system (2.2) does not undergo 
P-bifurcation. 

Proof. When JLL = 0, // 7 = 0. Then system (2.2) becomes 


dr = 

f /O 

r 

dt + 

( u 3 

t± r 2 


DC 



S 8 y 


dW r . 


(3.1) 


When ju 4 — 0 , equation (3.1) is a determinate system, and 
there is no bifurcation phenomenon.. Here we discuss the 
situation jU 4 ^ 0 , let 

l 

s / \ 

M _Ml 
8 16 


m(r) = 


M 


cr{r) = 




8 


r. 


The continuous random dynamic system generated by (3.1) is 

t i 

(/>(t)x = x + ^ m^^x^ds + Ja(y . 


where ° dW r is the differential in the meaning of 
Stratanovich, it is the unique strong solution of (3.1) with 
initial value X . And m — 0, cr(0) = 0 , so 0 is a fixed point 
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of (p . Since m{r) is bounded and for any r -p 0 , it 
satisfies the ellipticity condition: cr(r) =£ 0 ; it is assured that 
there is at most one stationary probability density. According 

to the It o equation of amplitude r(t ) , we obtain its FPK 
equation corresponding to (3.1) as follows 

dp _ d 


dt dr 


dp 


\\( .Mi') 


1 

L 52 1 

Ma 2 

l 

ILl 8 J 

r 

7 

dr 2 { 

— -r 
_ 8 

p ! 


(3.2) 


(3.3) 


Let — = 0 , then we obtain the solution of system (3.2) 
dt 

pit) = c Icr -1 (7)1 exp f | 

The above dynamical system (3.2) has two kinds of 
equilibrium state: fixed point and nonstationary Motion. The 
invariant measure of the former is S () and it’s probability 

density is 8 x . The invariant measure of the latter is V and 

it’s probability density is (3.3). In the following, we calculate 
the lyapunov exponent of the two invariant measures. 

Using the solution of linear Ito stochastic differential 
equation, we obtain the solution of system (3.1). 

r(f) = ,'(0) 

r\ i i i i /r i i i i ^ ^ 

/(o)+ 


exp 


cr 


m 


(°M°) 


els + ^<j' {Qt)clW r . 


VoL J o y 

(3.4) 

The Lyapunov exponent with regard to // of dynamic system 
(p is defined as: 


^(//) = lim-ln||r(f)|| 


t->w f 


(3.5) 


Substituting (3.4) into (3.5), note that cr(0) = 0, cr'(o) = 0 , 
we obtain the Lyapunov exponent of the fixed point: 

W) 

\( ' r A 

In ||r (0)|| + m (0) J ds + a' (0) J dW r (s) 


= lim 

(->+co f 


Wr(t ) 


(3.6) 


= m'(0) + cr'(0) lim 

= m’( 0) 

- „ , M^_M± 

8 16 

For the invariant measure which regard (3.4) as its density, 
we obtain the Lyapunov exponent: 

i 1 

/^(v) = lim -J[m'(r) + <T(r)<r'(r)]cfc 

oo £ 


-L 


m (r) + 


cr(r)cr ( r ) 


p ( r ) dr 


-2f 

JR 


m 


(dT, 


cr 


(r) 


p ( r) dr 


( \2 

16 t \ 

(rj exp 

— m\r ) 


Ma 


= -32sl2jU. 


3 f 
2 
A 


M\ + 


exp 


v 

16 ^ 


th_tL i 

8 16 


M* 




\ 


8 16 


(3.7) 


Lt 1 Li. 

Let a = LL -\ — . We can obtain that the invariant 

1 8 16 

measure of the fixed point is stable when a < 0 , but the 
invariant measure of the non-stationary motion is stable when 
a > 0, so a = a D = 0 is a point of D-bifurcation. 

Simplify Eq.(3.3), we can obtain 

2(8/;, +/2 2 -/j 4 ) 

p st (r) = cr A4 , (3.8) 

where c is a normalization constant, thus we have 


St 


(r) = o(r v ) r 0, 


(3.9) 


2(8 p,+p^-p.) , 

where v = — . Obviously when V < — 1 , 

Ma 

that is n + Tfl _ Tfl < 0 , p (r) is a 8 function, when 
1 8 16 


8 16 

point of p st (r) in the state space, thus the system undergoes 


— 1 < V < 0 , that is n + a-L — — 1 >0 r = 0 is a maximum 
8 16 


jU 2 JU, 

D-bifurcation when V = — 1 , that is Li, -I = 0 , is 

1 8 16 

the critical condition of D-bifurcation at the equilibrium point. 
When V > 0 , there is no point that makes p st {r) have 

maximum value, thus the system does not undergo 
P-bifurcation. 

Theorem 3.2(Stochastic Pitchfork bifurcation) When 
jU 3 —0,p 7 < 0 . Then the system (2.2) undergoes stochastic 
pitchfork bifurcation. 

Proof. When ju 3 = 0, // 7 ^ 0 . then Eq(2.11) can rewrite 
as follows 

l 


dr = 
Let (j) = 


Mi + 


V 


Mi_ 

8 


Ml 3 

r -l — -r 

8 


dt + 


M± r 2 
8 


dW r .(3.10) 


Mi 


r, p 7 <0 ,then the system (3.10) becomes 


f Vi) 

<p- 

-<p 3 

dt + 

f 

LI 8 J 




L 8 y 


d(p = 
which has the solution 


(p°dW, (3.11) 
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t/texp 


Mi 


Mi 


t + 


f/O 
U J 


W\(o) 


r 


1 + 2^ 2 J 


exp 


f 


2 

V V 


M\ 


Mi 

8 


t + 


M,' 

8 J 


W s (co) 


ds 


22 2 

(3.12) 

We now determine the domain D u (t,co) , where 
D fi (f , co) \cp e 91 : (t, co, cp) e Z)} (Z) = 9? x Q x X ) 

is the (in general possibly empty) set of initial values (j) e 91 
for which the trajectories still exist at time t and the range 
R u , (t, o>) of V Ml ( t,a > ) : D m (*,»)-> R p (t, co) . 

We have 

f 91, t> 0, 

D (t,co) = \ ( , s t s\ (3.i3) 

|(-J M (t,co),d Mi (t,co)), t< 0, 

where 

d f „ (*,<») = 


f 

t 

f 

i x 


X 

2 

V 

Jexp 

0 

2 i + i 

V 

y 

c/s 

y 




> 0 , 

and 

(- r A t ’ a, )’ r A t ’ a, ))’ t> °> 

9?, t < 0, 

where 

/ i 

r u , v r ,0 

exp 


(3.14) 


« + f Ht) V ' W 


exp 




Mi + 
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Z + 2 


Ma_ V 
v 8 J 






V 

> 0 . 

We can now determine 

:= ^.R (*> *>) 

and obtain 


M®)= 


( y, Ai + ^ > 0. 

{0}, 




where 

0 <d* ( t,CO ) = 




exp 
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2 

V V 
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Mi + 
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XA 


v 


t + 
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ds 


22 


V 

< oo 

The ergodic invariant measures of system (3.10) are 

fdj . /A 

(i) For //j H — < 0 , the only invariant measures is M m 

8 

=*<>■ 

(ii) For /j, + > 0 we have the three invariant forward 

8 

Markov measures jU^' = S Q and = 8 ±k ^ , where 


k *{°>)-= 


2 I 


exp 


x 


Mi + 


M h 

8 


t + 2 


M ± 

V 8 y 


'w;W 


X 


ds 

2 2 


We have Zi^ = OC. Solving the forward Fokker 
Planck equation 


L p = — 

I Z/i 




^(«* 2 AW) =0 


16 


yields 

(i) P M =^ 0 f °rallA +Y’ 

(ii) for p Ui > 0 

W = 


2(M+^) 
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exp 




0, 


^>0, 

^<0, 


and W = (- ^)> where ^ = T 
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...... _ e Fokker-Planck equation 

Naturally the invariant measures V ± — 0 + , are 


±*«H' 

those corresponding to the stationary measures q + v Hence all 

invariant measures are Markov measures. 

We determine all invariant measures(necessarily Dirac 
measure) of local RDS X generated by the SDE 


d(j) = 


A + 


A 




dt + 


on the state space 91 , u, + — - e 91 and 

8 

now calculate the Lyapunov exponent for each of these 
measure. 

The linearized RDS x t — DY ( t, 60, (j)) X satisfies the 
linearized SDE 


^ 4 v 


v * y 


f /^ 
V 8 j 


d%, = 


A + 


A 


-3(T (t, a), </>))' x, 


dt 




X, °dW. 


V8y 

Hence 

DY{t,co,(j))x = ^exp 


A + 


A 


f + 


f/O 
U J 


1 


l 

W t ds 


^(vr)=A 


8 


A 


So Vj 1 is stable for /u x H < 0 and unstable for //j + 


8 


— > o. 

8 


A 


A y 


L =- 


A + 


A 


vvv 




2 y 


A /A 

16 


(^pM) =0 ’ 


(zWW. (3.16) 


which has the unique probability density solution 


2 | rt+v- 


-1 


> 0 ■ We 


P Mi {(/)) = N fl t/) /,; exp 

Since 




, ^ > 0 . 


E^ q (p 2 = E(d^‘)“ =|(pV 1 ((p)d(p<oo, 

0 

the IC is satisfied. The calculation of the Lyapunov exponent 
is accomplished by observing that 

A ( u t coy = 


f 


exp 


f 


A + 


A 


t + 2 


^a¥ 
V 8 J 


A 




Y— invariant measure, its 


V 

Thus, if V r = 8 , 

Lyapunov exponent is 

A (/./) = lim-log | DY(t, co, 0)|| 

t 

= n x + — -31im|(Y(y,&>,^))" ds 
8 0 

= M +^-3E^ 0 2 , 

O 

provided the IC (f~ e L 1 ( P ) is satisfied. 

fi) For //j + £ 91 the IC for = 8 0 is trivially 

8 

satisfied and we obtain 

A 


2 J exp 

-00 

y(>) 

2V 

where 

/ 

v F(f)= Jexp 


A + 


A 


s + 2 


r LhX 2 

V 8 J 


w.H 


ds 


/in 

A + — 
1 8 


s + 


8 / 


W,(a) 


ds. 


Hence by the ergodic thoerem 

E ( rf - 1 ) 2 = ^ 7 lo s ^ ( ' ) = « + A • 


8 


finally 

A')=-2 


A + 


A 


A 

8 


< 0 . 


(iii) For //, H > 0, V 2ro = <5 M is F_ x measurable. 

8 ’ “ 


Since 


L 
thus 


(d^‘ ) = L (d* ) ,E (-d^ 1 ) 2 = E (d^ 1 ) 2 = p, 


p 2 


Pi + 


< 0 . 


^(v?) = -2 

v ^ y 

The two families of densities (r/ ;/ j >0 clearly undergo 


(ii) For //, H > 0, V 2 ' (y = <5 M is F _ rjj measurable, 

8 ’ d<0 

hence the density p Ml of p Ml = Ev 2 ‘ satisfies the 


Mi 

a P-bifurcation at the parameter value /l u , = Hence, we 
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have a D-bifurcation of the trivial reference measure 5 0 at 

/U D = 0 and a P-bifurcation of p = 

IV. P-BIFURCATION 



In the following, we consider the steady-state probability 
density p a ( r j of the linear It o stochastic differential equ- 
ation. Calculating extreme values of the invariant measure is 
one of the most popular efficient methods in studying the 
bifurcation of a nonlinear dynamical system. The invariant 
measure is an important characteristic value of stochastic 
bifurcation. 

4. 1 Case I: jU 3 ^ 0 , jU 7 =0 






When jU 3 ^ 0, // 7 = 0 . According 
amplitude r(t ) , we obtain its FPK equation as follows 

, Mi x Mi 

8 r 





(4.1) 


with the initial value Condition // 7 = 0, p{ /', / 1 r Q ,t Q ) - 

8 (r — /') j , / — > t 0 , where p(r,t\r Q ,t 0 ) is the transitio 
probability density of diffusion process r(f) ■ The invari; 
measure of r(f) is the steady-state probability density 
p st ( r) which is the solution of the degenerate system as 
following 

0 = - — 

dr 


(4.: 



Through calculation, we can obtain 

3 

Mi 



V 

f 

— V 

r(2-v) 

r 

y 

V 

J ) 


v—2 


(4.3) 

where V = (8/V, + // 2 ) // ' , T(x) = £ t x X e 'dt. 

According to Namachivaya’s theory[27], the extreme 
value of an invariant measure contains the most important 
essence of the nonlinear stochastic system. In other words, 
the invariant measure can uncover the characteristic infor- 
mation of the steady state. When the intensity of noise tends 
to zero, the extreme values of p st (r) approximately show 

the behavior of the deterministic system. If the process r (?) 


N. 


tf 


r.i|!irn 
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is ergodic then p t ( /■) can be regarded as the time 

measurement for staying in the neighborhood of aft) 

according to Oseled- 
ec ergodic theorem. 

From the analysis above, we know that the parameters fl 3 
> 0, /A > jU 4 >0. If p st \r) has a maximum value at r , 
the sample trajectory will stay for a longer time in the neighb- 
orhood of r , i.e. r is stable in the meaning of probability 
(with a bigger probability). If p si ( r ) has a minimum value 
(zero), it is just the opposite. 

* 

We now calculate the most possible amplitude r of 
system (2.9)., i.e. p t ( rj has a maximum value at r .So 
we have 


d P st ( r ) 


dr 


0, 


d 2 P st ( r ) 


dr 


,<0 


and the solution r = 0 or 


r = r = 


-8 Mi 


%M\ + Mi~ Ma 

Further, we have 

d 2 P s ,{r ) 


f 8 //, + 1 ^ 

as — — — < - 

Mi 2 


dr 


r= 0 " 


26+3(8 p x +p 2 -p 4 )p 1 ^2+{?,p i +p 2 -p A )p 1 > Q 


c £pJj) 

dr 


(8/f, + Mi M 4 ) 


8 M 3 


8 M 3 M 4 

8 Mi +Mi ~Mi 


Mi 


-16(8 p l +/u 2 -/j 4 f 


< 0 . 

• 4 s fV / \ % 

Thus what we need is r = r . In the meantime, p t tr ) is 0 
(minimum) at r = 0, . This means that the system subjected 
to random excitations is almost unsteady at the equilibrium 
point (r = 0) in the meaning of probability. The conclusion 
is to go all the way with what has been obtained by the 
singular boundary theory. The original nonlinear stochastic 
system has a stochastic Hopf bifurcation at r = r. 

2 , 2 — 8 Mi (■ 

x l +x 2 = £ , (i.e. r = r). 

8 M^Mi-Mi 

We now choose some values of the parameters in the equa- 
tions, draw the graphics of p _ t ( r j . The curves in the graph 

belonging to the condl,2,3,4 in turn are shown in Fig. la. It is 
worth putting forward that calculating the Hopf bifurcation 
with the parameters in the original system is necessary. If we 
now have values of the original parameters in system (2.1), 
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that b = 1. 5, ju = 1, co -0.8, r a =0.5625, x a =0.22, 
C, = C =0.1, y x =1.3,/c = 1.5,<5 A =0.1, B = l,x 0 = 
0. 1, Oj = 0. 1, (7 2 =0.1. After further calculations we 
obtain //, =-0.7625, // 2 =4.26167, //, =0.324941, // 4 
= 2.56999, 

v = 8//| = -0.715307 < - , 

Ma 2 

25.3905r 2 

P 1 ] ~ (2.59953 + 2.56999r 2 ) 271531 ' 

What is more is that r = 0.76791 1 where p st (r) has the 
maximum value(see Fig. lb). 

4.2 Case II: JU 3 — 0, /f 7 ^ 0. 

When jU, = 0, // 7 ^ 0. then Eq(2.9) can rewrite as 
ing 


follo\ 


dr = 


A 


A 

8 


, Ml 3 

r H — -r 

8 


dt + 


i 

/■ ,, 32 

— r 2 (4.4) 


V 


8 


J 


A 

According to the Ito equation of amplitude r(t'j , we 
obtain its FPK equation form (4.4) as following 


dp 

dt 


M\ 


Mi 


Mi, 3 

r —r-r 


2 Mi 


P(r ) 


aa a! 

2p n dr 


-(r 2 p(r)) 


(4.5 

with the initial value condition fp =0, p{r, t\ r 0 ,t 0 ) - 
S(r — r 0 ), t — > f 0 , where /?(r, 1 1 r (p / (| ) is the transitioi 
probability density of diffusion process r(Y). The invari; 
measure of r(/) is the steady- state probability density 
Pst( r ) i s the solution of the degenerate system a 

following 



f /C 

A + — 
1 8 


2a 


P( ? ') 


Through calculation, we can obtain 

Pj/) = 


■ (4. 


exp 


f 2 3 , 

r Mi / 
Ma J 




Mi 


f Mi\ 

V 8 ) 

Ma 


Mi 


Ml 


Mi 

V Ma ) 


Mi 


According to Namachivaya’s theory[27], the extreme 
value of an invariant measure contains the most important 
essence of the nonlinear stochastic system. In other words, 
the invariant measure can uncover the characteristic infor- 
mation of the steady state. When the intension of the noise 

tends to zero, the extreme values of p st ( r ) approximate to 

show the behavior of the deterministic system. If the process 
r(t) is ergodic then p t ( /- j can be regarded as the time 

mea- surement for staying in the neighborhood of a(t ) 
according to Oseledec ergodic theorem. 

From the analysis above, if p t ( r j has a maximum value 


at r , the sample trajectory will stay for a longer time in the 

neighborhood of r* , i.e. r is stable in the meaning of prob- 
ability (with a bigger probability). If p has a minimum 

value (zero), it is just the opposite. 

We now calculate the most possible amplitude r of 
system (4.4)., i.e. p st ( r ) has a maximum value atr* . So 
we have 


dp s M) 

dr r=r * 


d 2 p st (r) 


dr 2 


r=r* 


<0 


and the solution r — r 


AA+ 4 A 

8v 2 


. The probabilities and 


the positions of the Hopf bifurcation occurrence with differe- 
nt parameter are listed, and the corresponding results can be 
seen in Fig.2 as well. 

Since 


d 2 Pst (r) 


dr 2 


r=r 


3 ( fl + J )/<2 

4 + 5 


= 2 


Mi 


exp 
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( Y 

4 + — 
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Mi 

f \ 

Mi 


l 8 

v >^3 / 
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V Ma J 



< 0 . 

Thus what we need is r = r. The conclusion is to go 
all the way with what has been obtained by the singular boun- 


N. 


rS 


r.i|!irn 

Kcscarch Publlullon 


48 


www.ijntr.org 


Noise and Delay Induced on Dynamics of the Aeroelasticity Model 


dary theory. The original nonlinear stochastic system has a 
stochastic Hopf bifurcation at r = r. 


8 


Mi + 


Mi 


\ 


Mi + 4 /i 4 


2 . 2 
*1 +X 2 =■ 


8 Mi 


P s t (r) 


Fig. 2a a<- 



P Bt (r) Fig. 2b a=^- 



P st (r) Fig. 2c a<~ 



Fig.2. P-bifurcation of p a (r) at a p — .where 


a = Mi 


Mi 


CT~ — — 


8 p n 

We now choose some values of the parameters in the 
equations, draw the graphics of p st ( r j . It is worth putting 
forward that calculating the Hopf bifurcation with the para- 


meters in the original system is necessary. If we now have 
values of the original parameters in system (2.3), that 

b = 1.5,ju = l,a = l,r a = 0.5625, x a = 0.19,^ = 
C =0.3,7, = 1.4, k = 1.5, =0.2, 5 = 1.1, x 0 = 

0.5, cr, =0.3, C 2 =0.2. After further calculations we 
obtain //, = -0.2, /i 2 = 9.5907 ,ju 3 = 0 ,ju 4 = 5. 1 183, . 
// 7 = —20.4732, then 

p(r) = 0.28209e^ r2 . 


What is more is that r = 0.28209 where p st (r ) has the 
maximum value. 

4.3 Case III: // 3 ^ 0, // 7 ^ 0. 

When jLL ^ 0, /J 7 ^ 0. Similar above Case II discussions. 
Through calculation, we can obtain 


p sl (r) = cexp 


f 2 A 

r Mi 

V Ma / 


2 j ^M\ , M 2 8 / 6/6 

(&ju 3 +r 2 ju 4 ) "4 ft ft (48) 

where C is a normalization constant, according to Namachi- 

vaya’s theory, we calculate the most possible amplitude r" 
of system (2.9), we have 
r = r = 


1 -8 A -M 2 +M 4 - V ( 8 A +M 2 -M 4 ) 2 - 3 2 MiMi 

i 2 m 7 

In the meantime, p st {r ) is 0 (minimum) at r = O.This 
means that the system subjected to random excitations is 
almost unsteady at the equilibrium point ( r = 0 ) in the 
meaning of probability. The conclusion is to go all the way 
with what has been obtained by the singular boundary theory. 
The original nonlinear stochastic system has a stochastic 

Hopf bifurcation at r = r . 

1 , 2 

X l +X 2 = 

^ Mi Mi Mx \j ^ ^Mi Mi M \ ) 32//,// 

2 Mi 

(i.e. r = r ). 

We now choose some values of the parameters in the 
equations, draw the graphics of p sl (r) ■ It is worth putting 

forward that calculating the Hopf bifurcation with the 
parameters in the original system is necessary. If we now 
have values of the original parameters in system (2.6), that 

^ = ll,// = 1.2,® = l,r a =0.5,^ =0.18,^ =C = 

0.2, 7, = 1, k = 1 . 1, 8 a = 0.2, B = 1, x Q = 0.5, cr, = 0. 1, 
<J 2 = 0.3. After further calculations we obtain //, = —0.46 
25, // 2 =18.9095, >u 3 =0.510106, p 4 =12.284 ,// 7 =- 
0.3, 
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155.736e~ a227225r x 2 
P r _ (2.68 + 0.990208r 2 ) 4 ' 40395 ' 

What is more is that r = 0. 162777 where p a (r \ has the 
maximum value. 
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